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We have generalised the exact solution of the Riemann problem in special relativistic 
hydrodynamics (Marti & Miiller 1994) for arbitrary tangential flow velocities. The solu- 
tion is obtained by solving the jump conditions across shocks plus an ordinary differential 
equation arising from the self-similarity condition along rarefaction waves, in a similar 
way as in purely normal flow. The dependence of the solution on the tangential veloc- 
ities is analysed, and the impact of this result on the development of multidimensional 
relativistic hydrodynamic codes (of Godunov type) is discussed. 



1. Introduction 

The decay of a discontinuity separating two constant initial states (Riemann problem) 
has played a very important role in the development of numerical hydrodynamic codes 
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in classical (Newtonian) hydrodynamics after the pioneering work of Godunov (1959). 
Nowadays, most modern high-resolution shock-capturing methods (LeVeque 1992) are 
based on the exact or approximate solution of Ricmann problems between adjacent nu- 
merical cells and the development of efficient Riemann solvers has become a research 
field in numerical analysis in its own (see, e.g., the book of Toro 1997). 

Riemann solvers began to be introduced in numerical relativistic hydrodynamics at 
the beginning of the nineties (Marti, Ibahez & Miralles 1991). Presently, the use of high- 
resolution shock-capturing methods based on Riemann solvers is considered as the best 
strategy to solve the equations of relativistic hydrodynamics in nuclear physics (heavy 
ion collisions) and astrophysics (stellar core collapse, supernova explosions, extragalactic 
jets, gamma-ray bursts). This fact has caused a rapid development of Riemann solvers 
for both special and general relativistic hydrodynamics (see, e.g., the reviews by Ibahez 
& Marti 1999, Marti & Miiller 1999). 

In a previous paper (Marti & Miiller 1994; Paper I in what follows), we derived the 
analytical solution of the Ricmann problem for an ideal gas in special relativistic hydrody- 
namics for initial states where the flow is normal to the initial discontinuity. This solution 
has proven to be a useful for (i) the generation of analytical solutions to test relativistic 
hydrodynamic codes, and (ii) the development of numerical hydrodynamic codes based 
on an exact Riemann solver (e.g., Marti & Miiller 1996; Wen, Panaitescu & Laguna 1997). 
Because the solution only holds for flows which are normal to the initial discontinuity, 
numerical simulations based on the exact Riemann solver of Marti & Miiller (1994) are 
restricted to one-dimensional (ID) flows. Balsara (1994) and Dai & Woodward (1997) 
have circumvented this restriction at the price of constructing multi-dimensional Go- 
dunov schemes using a Riemann solver based on the two-shock approximation where the 
rarefaction wave is treated as a shock wave. An iterative relativistic nonlinear Ricmann 
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solver that takes into account the effects of nonvanishing tangential velocity components 

has already been implemented by Falle & Komissarov (1996) for Riemann problems in- 
volving two strong rarefactions. Finally, we note that the effect of the tangential velocity 
on the properties of simple waves and shocks has usually been ignored in relativistic 
hydrodynamics, because it is always possible to choose a reference frame in which the 
tangential velocity vanishes on both sides of the simple wave or shock. In the case of a 
Riemann problem this reference frame does not exist. The general solution cannot be 
constructed in terms of the solution of the purely normal flow case and by means of a 
single Lorentz transformation. 

In the following we derive the exact solution of a Riemann problem in Minkowski space- 
time with arbitrary tangential velocities. The solution can be implemented in multidi- 
mensional special relativistic hydrodynamic codes based on directional splitting, because 
it allows the computation of the numerical fluxes at every zone interface. Furthermore, 
according to recent work by Pons et al. (1998), who solve the equations of general rela- 
tivistic hydrodynamics using special relativistic Riemann solvers, this exact solution can 
also be implemented in general relativistic hydrodynamic codes. The exact solution also 
allows one to test the accuracy of other approximate Riemann solvers and codes. 

In this paper we closely follow the structure and notation used in Paper I, where 
one can also find the basic references for the theory of relativistic simple waves and 
shocks, first discussed by Taub (1948). Two key references in the theory of relativistic 
fluids are the review by Taub (1978) and the book by Anile (1989). The main idea 
behind the solution of a Riemann problem (defined by two constant initial states, L 
and R, left and right of their common contact surface) is that the self-similarity of the 
flow through rarefaction waves and the Rankinc-Hugoniot relations across shocks allow 
one to connect the intermediate states /* (/ = L, R) with their corresponding initial 
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states, /. The analytical solution of the Ricmann problem in classical hydrodynamics 
(see, e.g., Courant and Friedrichs 1948) rests on the fact that the normal velocity in the 
intermediate states, vf , can be written as a function of the pressure pi t in that state 
(and the flow conditions in state I). Thus, once pj„ is known, and all other unknown 
state quantities of i* can be calculated. In order to obtain the pressure pi t one uses the 
jump conditions across the contact discontinuity, which are given by 

PL* =Pr*{=P*) (1-1) 

«2*(p.) = «S.(p.)- (i-2) 



Equation (1.2) is an implicit algebraic equation in and can be solved by means of an 
iterative method. The function uf»(p*) is constructed by using the relations across the 
corresponding wave connecting the states / and I*. 

In the case of relativistic hydrodynamics the same procedure can be followed, the major 
difference with classical hydrodynamics stemming from the role of tangential velocities. 
While in the classical case the decay of the initial discontinuity does not depend on the 
tangential velocity (which is constant across shock waves and rarefactions), in relativistic 
calculations the components of the flow velocity are coupled through the presence of the 
Lorentz factor in the equations. In addition, the specific enthalpy also couples with the 
tangential velocities, which becomes important in the thermodynamically ultrarelativistic 
regime. 

The structure of the paper is the following. First, in §2 we present the equations of 
relativistic hydrodynamics for a perfect fluid in 3 spatial dimensions. In §3 and §4 we 
summarise the properties of the flow across rarefaction waves and shocks, respectively, 
and explain how to obtain the pressure and velocities behind the corresponding waves as a 
function of the state ahead the waves. In §5 we combine the results from the two previous 
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sections to solve the Riemann problem. Finally, in §6 we discuss the implementation of the 

solution in numerical relativistic hydrodynamics. Throughout the paper we will recover 

the corresponding purely normal flow expressions whenever it is of interest. 

2. The equations of relativistic hydrodynamics 

Let J M and (fi, v = 0, 1,2,3) be the components of the density current and the 
energy-momentum tensor of a perfect fluid, respectively 

jv = pu» (2.1) 
T^ v = phuV + p^ u , (2.2) 

where p denotes the proper rest-mass density, p the pressure, h = 1 + e+p/p the specific 
enthalpy, e is the specific internal energy and w M is the four- velocity of the fluid, satisfying 
the normalization condition 

= -1 (2.3) 
(throughout this paper, we will use the summation convention over repeated indices and 
units in which the speed of light is set to unity). In Cartesian coordinates, x M = (t, x, y, z), 
the Minkowski metric tensor rf ,v is given by 

rr=diag(-l, 1,1,1). (2.4) 

The evolution of a relativistic fluid is determined by the conservation equation of 
rest-mass (continuity equation) and energy-momentum 

J£ = 0, (2.5) 

T£" = (2.6) 
where iM stands for the partial derivative with respect to coordinate The above 
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system is closed by an equation of state (EOS) which we shall assume as given in the 
form p = p{p, e). 

The normalization condition of the velocity (p~a) leads to 



u" = W(l,v x ,v v ,v z ) (2.7) 

where W, the Lorentz factor, is 

W= {l-v 2 )- 1/2 (2.8) 

and 

v 2 = {v x f + {v y ) 2 + {v z ) 2 . (2.9) 

The equations of rclativistic hydrodynamics admit a conservative formulation which 
has been exploited in the last decade to implement high-resolution shock-capturing meth- 
ods. In Minkowski space time the equations in this formulation read 

U t +F^=0 (2.10) 

where U and F^(U) (i = 1, 2, 3) are, respectively, the vectors of conserved variables and 
fluxes 

U= (D,S\S 2 ,S 3 ,t) t (2.11) 

PW = (Dt^SV +p5 ll ,S 2 v l +pS 2l ,S 3 v l +pS 3l ,S l - Dv l ) T . (2.12) 

The conserved variables (the rest- mass density, D, the momentum density, S l , and the 
energy density r) are defined in terms of the primitive variables,(p,v l , e), according to 

D = pW 

S l ^phW 2 v i (2.13) 
t = phW 2 - p-D. 



System (2.10) is closed by means of an equation of state that we shall assume as given 
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in the form 

p=p(p,e). (2.14) 
The sound speed, c s , is then defined by 

(2.15) 



■ 2 dp 

where s is the specific entropy. 

In the case of an ideal gas with constant adiabatic exponent, 7, that we have considered 
in all the tests shown in this paper, the equation of state is simply 

p=( 7 -l)pe. (2.16) 

The hyperbolic character of the equations of relativistic hydrodynamics for causal 
equations of state and their eigenstructure are well known (e.g., Anile 1989). The com- 
plex dependence of the characteristic fields on the tangential velocity in arbitrary (i.e., 
noncomoving) frames, which is explicitly known (e.g., Donat et al. 1998), strongly affects 
the eigensystem and determines the properties of the waves. 

3. Relation between the normal flow velocity and pressure behind 
relativistic rarefaction waves 

Rarefaction waves are simple waves in which the pressure and the density of a fluid 
element decreases when crossing them. Choosing the surface of discontinuity to be normal 
to the x-axis, rarefaction waves would be self-similar solutions of the flow equations 
depending only the combination £ = x/t. Getting rid of all the terms with y and z 



derivatives in equations ( 2.1C ) and substituting the derivatives of x and t in terms of the 



derivatives of £, the system reads 

(v x -O^HpW 2 v*(v*-o+p}^+pW 2 vy(v*-o^+ P w\*(v*-o^ = o (3.1) 
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phW^v*-O^ + (l-v*O^=0 (3.2) 

pWV-^-^=0 (3.3) 
phW*(y°- 0^-^| = 0. (3.4) 
Equation ( |3.l|) comes from the equation of continuity. The remaining equations come 



from the momentum conservation. Equations (p.3|)-(3.4) reflect the expected result: If 
the tangential velocities are zero in the chosen state, no tangential flow will develop 
inside the rarefaction. Finally, instead of considering the conservation of energy, we use 
the conservation of entropy along fluid lines (following the same reasoning as in Paper I) 
which provides a relation between dp/d£, dp/d^ and dh/d^ 

dp 2 dp dh 

d~r hCs d~r p dt- (3 - 5) 



Non-trivial similarity solutions exist only if the determinant of system ( |3.1| )-( |3.5| ) van- 
ish. This leads to the condition 



V x (l - C|) ± - ^)[1 - gg _ („s)2 (1 _ C 2) 

? = 1 5"! > ( 3 - 6 

1 — v z d 



the plus and minus sign corresponding to rarefaction waves propagating to the left 7£«_ 
and right 7£_> , respectively. It is important to note that the two solutions for £ correspond 
to the maximum and minimum eigenvalues of the Jacobian matrix associated to the 
fluxes FwfU) (see, e.g., Donat et al. 1998), generalizing the result found for a vanishing 
tangential velocity in Paper I (equation (32)). 



After some manipulation, the system (3.1)-(|3.5|) can be reduced to just one ordinary 



differential equation (ODE) and two algebraic conditions 

phW 2 (v x -Qdv x + (l-(v x )dp = (3.7) 

hWv v = constant (3.8) 

hWv z = constant, (3.9) 
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with £ constrained by (3.6). From equations (3.S) and (3.9) it follows that v v /v z = con- 
stant, i.e., the tangential velocity does not change direction along rarefaction waves and 
it is only allowed to change its absolute value. Notice that, in a kinematical sense, the 



Newtonian limit (v l <C 1) leads to W = 1, but equations (3.8) and (3.9) do not reduce to 
the classical limit v v,z — constant, because the specific enthalpy still couples the tangen- 
tial velocities. Thus, even for slow or moderately relativistic flows (W « 1), the Riemann 
solution presented in this paper must be employed for thermodynamically relativistic 
situations (h > 1). The same result can be deduced from the Rankine-Hugoniot relations 
for shock waves (see next section). 

Using (|3.6D, the ODE (|3.7[) can be rewritten as 



dv x 



± 



dp phW 2 c s y/l + 5 (£ ± , V x , XT) 



(3.10) 



where v f = (v v ) 2 + (v z ) 2 is the absolute value of the tangential velocity and 

9(£± ,„.,„.) = < "' )2({i - 1) 



x\2 ' 



(3.11) 



The sign of £± corresponds to the sign chosen in (3.6). In the limit of zero tangential 



velocities, ti* = 0, the constants in (3.8), (3.9) are zero and the function g does not 
contribute. In this limit and in case of an ideal gas EOS one has 



W 2 dv x = ±—dp = ±—dp, 



(3.12) 



IP P 
recovering expression (30) in Paper I. 

Considering that in a Riemann problem the state ahead of the rarefaction wave is 



known, the integration of (3.1C) allows one to connect the states ahead (a) and behind 
(b) the rarefaction wave. Moreover, using ( |3.6| ), the EOS, and the following relation 
obtained from the constraint KWv 1 — constant, 



vl = h a W a vV 1 ( < )2 



h 2 + {KW a vif 



1/2 



(3.13) 
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the ODE can be integrated, the solution being only a function of p b - This can be stated 
in compact form as 

v x b =K^( Pb ). (3.14) 
Function lZ°^.(p) is shown in figure 1, for different values of the tangential velocity 
v in state a, the various branches of the curves corresponding to rarefaction waves 
propagating towards or away from a. Rarefaction waves move towards (away from) a, if 
the pressure inside the rarefaction is smaller (larger) than p a . In a Riemann problem the 
state a is ahead of the wave and only those branches corresponding to waves propagating 
towards a in figure 1 must be considered. Moreover, one can discriminate between waves 
propagating towards the left and right by taking into account that the initial left (right) 
state can only be reached by a wave propagating towards the left (right). The presence 
of a tangential velocity in the limiting state restricts the value of the normal velocity 
within the rarefaction wave to smaller values as it is clearly seen from figure 1 . Once the 
pressure in the post-wave state has been obtained, the corresponding tangential velocity 



follows from (3.13) (see lower panel of figure 1). 



4. Relation between post-shock flow velocities and pressure for 
relativistic shock waves. 

The Rankine-Hugoniot conditions relate the states on both sides of a shock and are 
based on the continuity of the mass flux and the energy-momentum flux across shocks. 
Their relativistic version was first obtained by Taub (1948) (see also Taub 1978 and 
Konigl 1980). 

If S is a hyper-surface in Minkowski space time across which p, and T^ v are dis- 
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continuous, the relativistic Rankinc-Hugoniot conditions are given by 

[pU^Tlf, = , 



11 



(4.1) 



where is the unit normal to £, and where we have used the notation 



[F] =F a - F b , 



(4.2) 



(4.3) 



F a and F b being the boundary values of F on the two sides of £. 

Considering S as normal to the a;-axis, the unitarity of n v allows one to write it as 



^ = 1^04,1,0,0), 



(4.4) 



where V s is interpreted as the coordinate velocity of the hyper-surface that defines the 
position of the shock wave and W s is the Lorentz factor of the shock, 

1 



Equation (4.1) allows one to introduce the invariant mass flux across the shock 



j = W s D a {V s - v x a ) = W s D b {V s - v x b ). 



(4.5) 



(4.6) 



According to our definition, j is positive for shocks propagating to the right. Note that our 
convention differs from that of both Landau & Lifshitz (1987) and Courant & Friedrichs 
(1948), but is the same as in Paper I. 

Next, the Rankine-Hugoniot conditions ( |4. l[ ) , (fh^) can be written in terms of the 
conserved quantities D, S J and r, and of the mass flux as follows 



3 



' 1 " 

D 



[P] 



W 8 
D 



= 0, 



(4.7) 

(4.8) 
(4.9) 
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~S Z 



D 



= 0, 



(4.10) 
(4.11) 



We note that in deriving equations (4.7)— (4.11) we have made use of the fact that the 
mass flux is nonzero across a shock. The conditions across a tangential discontinuity 



imply continuous pressure and normal velocity (by setting j = in equations (4.7), (4.8) 



and (4.11), and an arbitrary jump in the tangential velocity. 



Equations (4.9) and (4.1C) imply that the quantity hWv v,z is constant across a shock 
wave and, hence, that the orientation of the tangential velocity does not change. The 



latter result also holds for rarefaction waves (see §3). Equations ( ff .7| ), (4.8) and (4.11) 
are formally identical to the corresponding equations in Paper I (equations (47)- (49)) and 
can be manipulated in the same way to obtain v x as a function of p b , j and V s . Using 
the relation S x = (r + p + D)v x and after some algebra, one finds 



v x = h a W a v x a + 



W s {p b -p a ) 



, W„v x 

h a W a + { Pb - p a ) 1 



PaW a 



(4.12) 



This expression looks like that obtained for vanishing tangential velocity, but its presence 



is hidden within the Lorentz factor in state a, W a . From (J4.9|) and (4.10) expressions for 



and vl can be derived 



KW a vf 



1 - (vtf 



-1 1/2 



(4.13) 



M + (h a w a vvy_ 

The final step is to express j and V s as a function of the post-shock pressure. From 
the definition of the mass flux we obtain 



(4.14) 



where (V^) corresponds to shocks propagating to the right (left). 

The Taub adiabat (Thorne 1973), which relates (only) thermodynamic quantities on 
both sides of the shock, and the EOS can be used to derive the desired expressions. The 
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Taub adiabat, the relativistic version of the Hugoniot adiabat, is obtained by multiplying 



(4.2) first by (/ut M ) a and subsequently by (hu^b, and by summing up the resulting 



expressions. After some algebra one finds 

[^ 2 ]=(- + -W (4.15) 

V Pb Pa J 

In the general case, the above non-linear equation must be solved together with the EOS 
to obtain the post-shock enthalpy as a function of pb- In the case of ideal gas EOS with 
constant adiabatic index, the post-shock density pb can be easily eliminated and the Taub 
adiabat can be rewritten in the form (see Paper I) 

• 2 A, (7-l)(Pa-P&)\ (7 - l)(Pa -Pb) , . K{p a -p b ) 2 



fcg 1 + w ^» _ w -jvu. r», h 'W. _ ^ = , (4.16) 

V 7Pfc / IPb Pa 

which is a quadratic equation for the post-shock enthalpy hb as a function of p&. One of 
the two roots is always negative and must be discarded as a physical solution. 



Next multiplying (4.2) by and using the definition of the relativistic mass flux (4.6) 
one obtains 



which after using the EOS to eliminate pb and inserting the physical solution of ( 4.15 ) 
gives the square of the mass flux j 2 as a function of pb- 

Using the positive (negative) root of j 2 for shock waves propagating towards the right 



(left), equation ( 4.17 ) allows one to obtain the desired relation between the post-shock 



normal velocity v% and the post-shock pressure p&. In a compact way the relation reads 

v r b =S^(p b ). (4.18) 

The function SfL (p) is shown in figure 2 for several values of the tangential velocity 
u* in state a. Its various branches correspond to shock waves propagating towards or 
away from a. In order to select the relevant branch of the function S^(p) (figure 2) the 
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same argumentation as in the case of rarefaction waves can be used (see §3). When v% is 
known, (4.13) can be used to determine v y and v z in the post-shock state. 



5. The solution of the Riemann problem with arbitrary tangential 
velocities. 

The decay of an initial discontinuity gives rise, in general, to three elementary nonlinear 
waves (see, e.g., Landau & Lifshitz 1987). Two of them can be shocks or rarefaction waves, 
one moving towards the initial left state and the other towards the initial right state. 
Between them, two new states appear, namely L» and separated from each other 
through the third wave, which is a contact discontinuity moving along with the fluid. 
Across the contact discontinuity pressure and normal velocity are constant, while the 
density and the tangential velocity exhibits a jump. Accordingly, the time evolution of a 
Riemann problem can be represented as: 

I -> LW^L*CR*W^R (5.1) 

where W and C denote a simple wave (shock or rarefaction) and a contact discontinuity, 
respectively. The arrows («— / — ►) indicate the direction (left / right) from which fluid 
elements enter the corresponding wave. 

As in the Newtonian case, the compressive character of shock waves (density and pres- 
sure rise across the shock) allows us to discriminate between shocks {S) and rarefaction 
waves (1Z): 

{^*- (->) ) Pb<Pa 
(5.2) 
<S<- (_) , Pb> Pa 

where p is the pressure and subscripts a and b denote quantities ahead and behind the 
wave. For the Riemann problem a = L(R) and b = L*(iZ*) for and W— y , respectively. 
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Thus, the possible types of decay of an initial discontinuity can be reduced to: 

(a) I — > L <S<_ C R* R p L < Pl* = Pr* > Pr (5.3) 

(b) I —> L L* C R* TZ^ R p L < pl* = Pr* < Pr (5.4) 

(c) I — > L TZ^- C R* R p L > pl* = Pr* < Pr- (5.5) 
The solution of the Riemann problem consists in finding the intermediate states, 

and as well as the positions of the waves separating the four states (which only 
depend on L, L*, R„ and R). The functions W_ > and W*_ allow one to determine the 
functions v^{p) and respectively. The pressure p* and the flow velocity in the 

intermediate states are then given by the condition 

«£*(?)= «L(P)=<- (5-6) 

When and have been obtained the remaining quantities can be computed. The 
EOS gives the specific internal energy and the remaining state variables of the interme- 
diate state 7* can be calculated using the relations between J# and the respective initial 
state I given through the corresponding wave. In particular, the tangential velocities can 



be calculated from ( |3.8|) - (3J5) for rarefaction waves and from the Rankinc-Hugoniot 



jump conditions ( |4.9| )-( |4.10| ) for shock waves. Notice that the solution of the Riemann 
problem depends on the modulus of v l , but not on the direction of the tangential veloc- 
ity. Figure 3 shows the solution of a particular Riemann problem (Sod 1978) for different 
values of the tangential velocity v v = 0, 0.5, 0.9, 0.99. The crossing point of any two lines 
in the upper panel gives the pressure and the normal velocity in the intermediate states. 
The range of possible solutions in the (p, i^-plane is marked by the shaded region. 
Whereas the pressure in the intermediate state can take any value between pl and pr, 
the normal flow velocity can be arbitrarily close to zero in the case of an extremely rel- 
ativistic tangential flow. The values of the tangential velocity in the states L* and i?* 
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are obtained from the value of the corresponding functions at v x in the lower panel of 
figure 3. 

To study the influence of tangential velocities on the solution a Riemann problem, we 
have calculated the solution of a standard test involving the propagation of a relativistic 
blast wave produced by a large jump in the initial pressure distribution (see Marti & 
Muller 1999) for different combinations of v L and v R . The initial data are pl = 10 3 , 
PL = 1, v x L = 0; pr = 10~2, pa = 1, v R — and the 9 possible combinations of 
v L R = 0, 0.9, 0.99. The results are given in figure 4 and table 1. 

The speeds of the waves propagating to the left and right, respectively, tend to zero in 



the limit of v L R — ► 1. This result is easily deduced from (3.6) which gives the velocities 



of the head (£h) and tail (£t) of the rarefaction wave, and from (4.14) which gives the 
shock speed (V s ). These three velocities tend to the normal velocity ahead of the wave 



(v x ) when the corresponding tangential velocity (n*) tends to w^ oa . = \J\ — (v x ) 2 - 

Although the structure of the solution remains the same (left propagating rarefaction 
wave, right propagating blast wave) the values in the constant intermediate states change 
by a large amount. The results show (see figure 4 and table 1) that ^ remains constant as 
long as v\ is constant. The value of £/, only depends on the left state, and decreases with 
increasing v L . The pressure in the intermediate state, p», increases with increasing v R 
due to the larger effective inertia of the right state (because of the increase of the Lorentz 
factor of the right state) . The velocity of the shock is determined by the competing effects 
of an increased p» and the larger inertia of the right state. The shape of the function 
lZ1±(p) (see figure 1) and its dependence on v L shows that both and v x decrease with 
increasing v L . On the other hand, as a change in v L does not alter the thermodynamic 
state, the flow evolution accross the rarefaction wave is along the same adiabat for any 
value of v L . This explains why pl* gets smaller with increasing v L (consistently with the 
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pL* 



PR* 



P* 



V* 



0.00 0.00 
0.00 0.90 
0.00 0.99 



9.16e-2 1.04e+l 1.86e+l 0.960 0.987 -0.816 +0.668 
1.51e-l 1.46e+l 4.28e+l 0.913 0.973 -0.816 +0.379 
2.89e-l 4.36e+l 1.27e+2 0.767 0.927 -0.816 -0.132 



0.90 0.00 
0.90 0.90 
0.90 0.99 



5.83e-3 3.44e+0 1.89e-l 0.328 0.452 -0.525 +0.308 
1.49e-2 4.46e+0 9.04e-l 0.319 0.445 -0.525 +0.282 
5.72e-2 7.83e+0 8.48e+0 0.292 0.484 -0.525 +0.197 



0.99 0.00 
0.99 0.90 
0.99 0.99 



1.99e-3 1.91e+0 3.16e-2 0.099 0.208 -0.196 +0.096 
3.80e-3 2.90e+0 9.27o-2 0.098 0.153 -0.196 +0.094 
1.29e-2 4.29c+0 7.06e-l 0.095 0.140 -0.196 +0.085 



Table 1. Solution of the relativistic Riemann problem at t = 0.4 with initial data pl = 10 3 , 
PL = 1.0, v£ = 0.0, pr — 10~ 2 , pr — 1.0 and — 0.0 for 9 different combinations of tangential 
velocities in the left (v L ) and right (vr) initial state. An ideal EOS with 7 = 5/3 was assumed. 
The various quantities in the table are: the density in the intermediate state left (pl*) and right 
{pR*) of the contact discontinuity, the pressure in the intermediate state (p*), the flow speed 
in the intermediate state (vi?), the speed of the shock wave (V s ), and the velocities of the head 
and tail (£ t ) of the rarefaction wave. 

decrease of pressure). Thus, the net effect of tangential velocities on rarefaction waves is 
to evolve further along the adiabat, and to reach smaller intermediate pressure values. 

6. Conclusions 

We have obtained the exact solution of the Riemann problem in special relativistic 
hydrodynamics with arbitrary tangential velocities. Unlike in Newtonian hydrodynam- 



18 J. A. Pons, J. Marti and E. Miiller 

ics, tangential velocities are coupled with the rest of variables through the Lorentz fac- 
tor, present in all terms in all equations. It strongly affects the solution, especially for 
ultra-relativistic tangential flows. The specific enthalpy acts as a coupling factor, too. 
It modifies the solution in thermodynamically relativistic situations (energy density and 
pressure comparable to or larger than the proper rest-mass density) even in slow flows. 

Our solution has interesting practical applications. First, it can be used to check the 
different approximate relativistic Ricmann solvers developed by various authors in the 
last decade, and to test multi-dimensional hydrodynamic codes based on directional 
splitting. Second, it can be used to construct multi-dimensional relativistic Godunov 
type methods. The latter project is currently in progress and will be reported elsewhere. 
Finally, using the procedure described by Pons et al. (1998) the exact solution can be 
used as a building block in a general relativistic hydrodynamic code. 

The computational cost of the exact Ricmann solver derived above is comparable to 
the one presented by Marti & Miiller (1994) which is valid for purely normal flows. 
Hence, when our exact Riemann solver is applied to multi-dimensional flow problems, 
the difference in efficiency with respect to most linearized Riemann solvers is reduced. 

This work has been supported by the Spanish DGICYT grant PB97-1432, and by the 
agreement between the Spanish CSIC and the Max-Planck-Gesellschaft. J. A. P. acknowl- 
edges financial support from the Spanish Ministerio de Education y Ciencia (FPI grant). 
We acknowledge an unknown referee B for his interesting comments which helped us to 
improve the quality of the manuscript. 
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Figure captions 

Figure 1: Loci of states which can be connected with a given state a by means of rela- 
tivistic rarefaction waves propagating to the left (7£_) and to the right (7£_) and moving 
towards or away from a. Upper panel: pressure versus v x ; lower panel: tangential velocity 
versus v x . Solutions for different values of the tangential velocity v* — 0,0.5,0.9,0.953 
correspond to solid, dashed, dashed-dotted and dotted lines, respectively. The state a is 
characterized by p a — 0.6, p a = 1.0, and v x — —0.3. An ideal gas EOS with 7 = 5/3 was 
assumed. 

Figure 2: Loci of states which can be connected with a given state a by means of 
relativistic shock waves propagating to the left (<S^) and to the right (<S_>) and moving 
towards or away from a. Upper panel: pressure versus v x ; lower panel: tangential velocity 
versus v x . Solutions for different values of the tangential velocity v* — 0,0.5,0.8,0.865 
correspond to solid, dashed, dashed-dotted and dotted lines, respectively. The state a is 
characterized by p a = 0.25, p a = 1.0, and v x — —0.5. An ideal gas EOS with 7 = 5/3 
was assumed. 

Figure 3: Graphical solution in the (p, v x )-plane (upper panel) and in the (v t ,v x )- 
plane (lower panel) of the relativistic Riemann problem with initial data pl = 1.0, 
pL = 1.0, i>£ = 0.0; pr = 0.1, pr — 0.125 and v R = 0.0 for different values of the 
tangential velocity v* = 0, 0.5, 0.9, 0.999, represented by solid, dashed, dashed-dotted 
and dotted lines, respectively. An ideal gas EOS with 7 = 1.4 was assumed. The crossing 
point of any two lines in the upper panel gives the pressure and the normal velocity in 
the intermediate states. The value of the tangential velocity in the states L* and i?» is 
obtained from the value of the corresponding functions at v x in the lower panel, and 
Io gives the solution for vanishing tangential velocity. The range of possible solutions is 
given by the shaded region in the upper panel. 
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Figure 4: Analytical pressure, density and flow velocity profiles at t = 0.4 for the 

relativistic Riemann problem with initial data pl = 10 3 , pl = 1-0, w£ = 0.0; pr = 10~ 2 , 

pn = 1.0 and vfj = 0.0, varying the values of the tangential velocities. From left to right, 

?4 = 0,0.9,0.99 and from top to bottom v\ = 0,0.9,0.99. An ideal EOS with 7 = 5/3 

was assumed. 
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